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SOME FEATURES OF ARTIFICIALLY THICKENED FULLY DEVELOPED TURBULENT 
BOUNDARY LAYERS WITH ZERO PRESSURE GRADIENT 1 

By P. S. Klebanoff and Z. W. Diehl 


SUMMARY 

An account is given herein of an investigation conducted to 
determine the feasibility of artificially thickening a turbulent 
boundary layer on a flat 'plate. A description is given of sev- 
eral methods used to thicken artificially the boundary layer. 
It is shown that it is possible to do substantial thickening and 
obtain a fully developed turbulent boundary layer , which is 
free from any distortions introduced by the thickening process , 
and , as such , is a suitable medium for fundamental research. 
Measurements of mean velocity , spectrum of u-fluctuation, and 
intensity of u-fluctuation which served as criteria for determin- 
ing the state of the layer are presented. Some features of the 
fully developed turbulent boundary layer with zero pressure 
gradient are discussed. The mean flow in the boundary layer 
is compared with the logarithmic laws of mean-velocity distri- 
bution derived for pipes and channels , and the existing theories 
dealing with the shape of the spectrum in isotropic turbulence 
are applied to the spectrum measured in the outer portion of 
the boundary layer. It is also shown that the turbulence in 
the outer portion of the boundary layer has an intermittent 
character similar to that found for wakes and jets. 

INTRODUCTION 

The lack of a satisfactory theory for turbulent shear flow 
and the development of the hot-wire anemometer as a power- 
ful tool for turbulence research have directed attention to 
the quantitative measurement of the characteristics of tur- 
bulence in shear flow as a means of obtaining information 
on which to base a sound theory. The recent concept of 
local isotropy proposed by Kolmogoroff (references 1 and 2), 
that is, that, whatever the nature of the large-scale motions, 
the small-scale turbulence tends to be isotropic, has intensi- 
fied the experimental program with the hope that the exist- 
ing theories of isotropic turbulence may be useful. While a 
considerable amount of work has been done in various flow 
fields (references 3, 4, 5, and 6), attention is here centered on 
the boundary layer, not only because of its practical impor- 
tance but to provide information which can be compared 
with that obtained in other types of shear flow and thus 
show differences and likenesses and in particular show 
whether certain phenomena are universal. It is felt that 
the results of the measurements can best be compared with 
the theory of isotropic turbulence if such complicating fac- 
tors as pressure gradient, curvature, surface roughness, and 
compressibility are eliminated. Furthermore, the turbu- 
lence should be fully developed and the Reynolds number 
should be sufficiently high to bring the domain of local iso- 


tropy within the range of measurement. In addition the 
boundary layer should be as thick as possible in order to 
facilitate the use of multiwire hot-wire probes and minimize 
errors due to wire length. 

These requirements mean that the surface should be 
smooth and flat and that high Reynolds number should be 
obtained by thickness rather than by high speed. Such 
requirements are difficult to meet, except in wind tunnels 
especially constructed for the purpose. The use of one wall 
of a conventional wind tunnel does not appear to be a sat- 
isfactory solution because of the likelihood of extraneous 
effects. The problem of obtaining a thick boundary layer 
was encountered at the National Bureau of Standards when 
a turbulent-boundary-layer investigation using a flat plate 
was undertaken. The test section of the tunnel was 19 feet 
long, but it was not feasible to use a plate longer than 12. 
feet. It was felt that the need for a thicker boundary 
layer than one which would normally develop on 12 feet of 
surface was sufficiently urgent here, and possibly in other 
laboratories with even less length available, to warrant an 
investigation of the possibilities of artificially thickening the 
layer. This involved not only methods of thickening but 
also setting up criteria by which to establish the identity of 
a fully developed layer over a smooth surface. 

In determining the feasibility and validity of artificial 
thickening a number of measurements had to be made, 
including mean velocities, intensity of w-fluctuations, and 
the spectrum of ^-fluctuations. As a consequence it is pos- 
sible to show how such quantities are affected by thickening 
devices and what values they assume when the layer is no 
longer influenced by its beginnings. It is obviously very 
necessary to know whether there exists a “flat-plate boundary 
layer” of universal character depending only on Reynolds 
number. This is shown to be the case insofar as the observed 
quantities provide a proper sampling. 

The investigation was conducted under the sponsorship 
and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

The authors wish to express their appreciation for the con- 
stant advice and active interest of Dr. G. B. Schubauer both 
in the experimental program and in the preparation of the 
report. Acknowledgment is made to Messrs. G. H. Adams 
and L. M. Kelly for their assistance in the experimental 
program. The authors also wish to thank Dr. C. M. Tchen 
and Dr. J. Laufer for their interest and stimulating discus- 
sions. 


1 Supersedes NACA TN 2475, "Some Features of Artificially Thickened Fully Developed Turbulent Boundary Layers with Zero Pressure Gradient" by P. S. Klebanoff and Z. W. 
Diehl, 1951. 
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distance along surface from leading edge of 
flat plate 

distance normal to surface, measured from 
surface of flat plate 
mean velocity in boundary layer 
mean velocity in free stream, just outside 
boundary layer 

x, y, and z components of instantaneous 
turbulent-velocity fluctuations 
mean-square values of u, v, and w 
root-mean-square value of u 
density of air 
kinematic viscosity of air 
static pressure 

static pressure at 2=10.5 feet 
free-stream dynamic pressure 
free-stream dynamic pressure at 2 = 10.5 feet 
dynamic pressure in boundary layer as 
measured with a surface tube 
shearing stress at surface 
absolute constant 
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longitudinal scale of ^/-fluctuation 
microscale in isotropic turbulence 
total dissipation of turbulent energy per 
unit mass per second 
time 

frequency, cycles per second 
wave number in intermediate range where 
initial forces and viscosity both play role 


one-dimensional wave number 



three-dimensional wave number 
normalized one-dimensional spectral func- 
tion associated with frequency n 


z s constant 

F(z) = UF(ri)/L x 

E x (k x ) one-dimensional spectral function associated 

with k x 

E(k) three-dimensional spectral function asso- 

ciated with k 


APPARATUS AND PROCEDURE 

GENERAL ARRANGEMENT 

The investigation was conducted in the 4^-foot wind tunnel 
at the National Bureau of Standards, using a flat plate 12 
feet long mounted vertically along the center line of the test 
section, the leading edge being 6 feet from the upstream end 
of the test section. The plate was ^-inch aluminum with a 
smooth surface and a sharp leading edge formed by a smooth 
symmetrical taper beginning 4 inches from the leading edge. 
The plate was fastened along the horizontal edges to one side 
of 3-inch channel irons anchored to the upper and lower tun- 
nel walls, and so divided the test section into nearly equal 
passages. It was considered desirable to conduct the inves- 
tigation in a tunnel of low turbulence so that the turbulence 
of the free stream would be negligible compared with that in 
the boundary layers. The turbulence level of the tunnel 
was 0.02 percent at 30 feet per second and 0.04 percent at 100 
feet per second. An elevation view of the tunnel is shown in 
figure 1. More details as to the arrangement of damping 
screens used in obtaining the low turbulence are given in 
reference 7. 



Figure 1.— Elevation view of NBS 4^-foot wind tunnel. 
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In order to obtain a condition of zero pressure gradient, a 
flexible second wall of sheet aluminum was mounted on the 
tunnel wall facing the working side of the plate. It was 
attached to the ends of screws threaded into the tunnel wall 
and thus provided an adjustable section to make the passage 
between the plate and the wall sufficiently diverging to offset 
the pressure drop resulting from boundary-layer growth. 

Another important consideration was proper angle of 
attack to afford smooth flow' entry at the leading edge on the 
working side of the plate. This was accomplished by dis- 
placing the upstream edge of the plate and by partial block- 
ing of the passage on the working side. Since the obstruc- 
tions causing the blocking had to be placed farther down- 
stream than any usable portion of the plate, the passage was 
extended 4 feet by joining the trailing edge of the plate with 
a section of plywood. 

Under these conditions the pressure distribution along the 
plate was found to be as shown in figure 2(a). From 2 feet 
on downstream the pressure is seen to be uniform on the 
average with variations of about 1 percent of q 0 about the 
mean. These variations could not be removed by the 
adjustable wall because they were associated with a slight 
waviness of the surface of the plate. 

Figure 2(b) shows the transition region at a speed of 108 
feet per second for smooth entry conditions at the leading 
edge. The beginning of transition corresponds to a Reynolds 
number (based on length of laminar flow) of 3,175,000, and 
here transition results from the growth of more or less regular 
laminar-boundary-layer oscillations. 



0 2 4 6 8 10 12 


ft 

(a) Pressure distribution along flat plate, 

(b) Variation of dynamic pressure as measured with a surface tube through transition 
region. Ui = 108 feet per second. 

Figure 2.— Pressure-distribution and transition data. 


INSTRUMENTATION AND TECHNIQUES 

Traversing mechanism. — In order to traverse in a hori- 
zontal plane normal to the surface the unit shown in figure 
3 was used. A micrometer screw with 0.5-millimeter pitch 
is driven through a set of beveled gears and a shaft by a 
selsyn motor remotely controlled. To provide for ease of 
alinement, the pressure probes and hot-wire holders are 
attached so that they can move fore and aft and rotate in 
a vertical plane. By means of the bracket shown in the 
figure the unit is bolted to a streamlined strut so that the 
end opposite the gears is forced against the surface. Channel 
irons anchored to the floor and ceiling of the tunnel served 
as guide rails for the strut, and traversing longitudinally was 
done by moving the strut fore and aft by means of a chain 
and sprocket manually operated from outside the tunnel. 

Measurement of pressure distribution, transition, and 
mean velocity. — Measurements of pressure distribution and 
transition were made with the pitot-static head shown in 
figure 4. The impact tube makes contact with the surface, 
and the static tube is opposite the point of contact inch 
away from the surface. The static and impact tubes are 
made of 0.04-inch nickel tubing with 0.003-inch wall thick- 
ness. The static tube is made according to the conven- 
tional design for such a tube, and the impact tube is flattened 
at the end to form a rectangular opening 0.014 inch wide. 
The static tube was used for the measurement of pressure 
distribution and both impact and static tubes for the meas- 
urement of transition. The head was attached by a long arm 
(containing the pressure leads) to the strut and could be 
moved back and forth in the ^-direction while the forklike 
guides maintained proper relation to the surface. The deter- 
mination of the transition point with the surface tube depends 
on the increase in impact pressure close to the surface when 
the layer becomes turbulent. The abruptness of the increase 
is shown in figure 2(b). The beginning and ending of the 
transition region coincide with the beginning and ending of 
the rise. This was checked with the hot-wire head shown 
in figure 4. This sledlike head was attached to the strut so 
that it could be moved fore and aft like the pitot-static head. 
The hot-wire w'as at a fixed distance of 0.03 inch from the 
surface, and visual observation of the fluctuations on an 
oscilloscope for different positions along the surface confirmed 
the indications of the surface tube. Measurements of mean 
velocity were made with a pitot-static combination similar 



Figure 3.— Traversing mechanism with total-head and static tubes attached. Reference 

scale in inches. 
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Figure 4— Sliding hot-wire head and surface pitot-static tube. 


to the one described above but without the guides. It is 
shown with the traversing unit in figure 3. No correction 
was made for the effect of turbulence on the readings. 

Measurement of intensity and spectrum of //-fluctuation. — 
A description of the method and the hot-wire equipment used 
for measuring the intensity of w -fluctuation is given in detail 
in reference 8. Platinum wire 0.0001 inch in diameter and 
0.5 millimeter long was used. This was wire drawn by the 
Wollaston process, and the method of making a hot-wire 
probe was to etch away the silver coating and then soft- 
solder the wire to the tips of fine sewing needles. The time 
constant of the wire was approximately 0.25 millisecond and 
was determined by the square-wave method described by 
Kovasznay (reference 9). Compensation for the time lag 
was made by a resistance-capacitance network in the ampli- 
fier. The response of the uncompensated amplifier was flat 
from 10 to 5000 cycles per second and fell off to about 
50 percent at 10,000 cycles per second. The mean-square 
output was read on a microammeter working out of a 
thermocouple circuit. 

In order to measure the spectrum of w -fluctuation the 
output circuit of the hot-wire amplifier was bypassed and 


the signal fed into an Electrical Research Products, Inc., 
sound analyzer with a frequency range from 10 to 10,000 
cycles per second and a choice of fixed band widths of nom- 
inally 5 or 50 cycles per second. The rectifier type of output 
circuit of the analyzer was replaced by an additional stage 
of amplification and a thermocouple and meter. The pur- 
pose' of this was to read mean-square currents and also to 
obtain a highly damped indicator. In reading the output, 
averages were taken over a period of 1 minute. The mean- 
square voltage corresponding to the compensated hot-wire 
signal passed by the band was obtained by first reading the 
hot-wire signal on the output meter and then matching this 
signal with a known sinusoidal signal fed into the system 
from a General Radio microvolter. The known signal had 
the frequency corresponding to the center of the band. The 
mean-square hot-wire voltage per unit frequency was then 
obtained by dividing by the effective band width. The 
effective band width was defined as the rectangular band 
width having the same area as the experimentally determined 
band shape. These were 

Effective Nominal 

(cps) (cps) 

7. 07 5 

57. 9 50 

It was found to be advantageous to use the 7.07 band from 
10 to 1000 cycles per second and the 57.9 band from 1000 
to 10,000 cycles per second. 

ARTIFICIAL THICKENING OF BOUNDARY LAYER 

It is at once obvious that objects protruding outward from 
a surface create more drag than the surface that they cover 
and hence offer a means of producing a layer of de-energized 
flow in a comparatively ^hort distance in the streamwise 
direction. The energy of the main flow is diverted into eddy 
motion in wakes, and, while the How is turbulent, it does not 
have the same dynamic structure as a fully developed 
turbulent boundary layer over a smooth surface. The ques- 
tion is, after leaving drag-producing objects, what length of 
smooth surface is needed for the layer to become the same 
as a smooth-surface generated layer. It would seem that 
the length would be less than that required to generate a 
layer of the same thickness, but this is by no means a foregone 
conclusion. 

Sinc£ the ultimate aim of the work was to produce a layer 
approximately 3 inches thick, this thickness having been 
chosen as sufficient for the convenient use of hot-wire probes, 
no effort was made to find the maximum thickness that could 
be produced in a given length. Emphasis was placed rather 
on simplicity of method and adequacy of tests. 

The turbulent boundary layer along the plate was first 
obtained by the conventional method of free transition, that 
is, from the instability of laminar-boundary-layer oscillations. 
Then, various devices such as rods, screens, and distributed 
roughness were used on the laminar layer to force transition 
and de-energize the flow. Attention was given to maintain- 
ing the condition of zero pressure gradient for the modified 
layers by adjusting the “false” wall. The first attempt at 
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thickening the layer was to use rods, which were 0.04, 0.11, 
0.18, and 0.25 inch in diameter. These were placed singly in 
the laminar layer and in contact with the surface at 4 feet 
from the leading edge. As will be seen, these were not satis- 
factory because of the length required for their distorting 
effect to disappear. The second method consisted in placing 
wire screens of various mesh sizes at different angles to the 
plate and close to the leading edge, so that the screen and 
plate formed a narrow wedge opening into the wind. A 24- 
mesh screen of 0. 0075-inch wire diameter, which was 1 foot 
wide, making an angle of 0.5° with the plate and extending 
from floor to ceiling, was a satisfactory thickening device, 
but it was not considered practical because of the possibility 
of accidental shifting. The third and final method consisted 
in cementing sandpaper to the surface beginning at the lead- 
ing edge and extending for distances up to 2 feet downstream. 
Various amounts of sand roughness were tried until the layer 
was of the desired thickness. The final arrangement of the 
roughness was No. 16 floor-sanding paper, cemented to the 
plate from the leading edge to 2 feet downstream and 
extending from floor to ceiling. 

Initially the one-dimensional spectrum of ^-fluctuations 
was thought to be the best over-all test to be applied to the 
resulting layers. Accordingly, for a basis of comparison, 
considerable attention was given to the spectrum of the fully 
developed layer resulting from free transition. It was found 
later that the mean-velocity profile was a still more sensitive 
test of the condition of the layer. Finally, the distribution 
of u f across the layer was added to the list of criteria. These 
three, compared in different parts of the layers and compared 
with other available data, served as the sample phenomena 
to determine the state of the layer. 

RESULTS AND DISCUSSION 

MEAN-VELOCITY MEASUREMENTS 

With free transition beginning at the 5-foot position, as 
shown in figure 2, mean-velocity distributions were deter- 
mined at 2=7.5, 8.5, 9.5, and 10.5 feet at a free-stream speed 
of 108 feet per second. These are shown in figure 5. A lower 
value of Ux would have shifted the transition point too close 
to the trailing edge of the plate. Figure 6 shows the same 
velocity profiles on a nondimensional basis. In order to avoid 
the arbitrariness of 5, the nondimensional distance from the 
wall y/8* is used. It is seen from figure 6 that downstream 
from the extreme position of the transition point, that is, 
2 = 0.75 feet, the profiles exhibit similarity. The two- 
dimensional nature of the flow was checked by making veloc • 
ity measurements 10 inches below and above the center line 
at 2 = 10.5 feet and comparing with the velocity distribution 
at the center line. Agreement among the three sets of 
distributions confirmed the two -dimensionality of the flow. 

The velocity measurements for the 0.04- and 0.25-inch 
rods are shown in figures 7 and 8. At a free-stream speed 
?7i = 108 feet per second, the flow, after passing over the 
rods, became turbulent in the separated layer and reattached 
to the surface in the turbulent state. At a lower free-stream 
speed U x = 30 feet per second for the 0.04-inch rod, visual 
observation of an oscilloscope trace obtained by traversing 


downstream from the rod at a distance of 0.03 inch from the 
surface, with a -hot-wire anemometer, showed the reappear- 
ance of laminar-boundary-layer oscillations, indicating that 
the separated layer returned to the surface as a laminar layer 
without the occurrence of transition. This effect was also 
noticed by Liepmann and Fila (reference 10). Figures 9 
and 10 give some of the representative profiles for the 0.04- 
and 0.25-inch rod in nondimensional form. The profile for 
2=5.17 feet in figure 9 differs from those at 6.0 . and 10.5 
feet by more than the experimental scatter. This would 
place the beginning of similarity for the 0.04-inch rod 
between 2=5.17 feet and 2 = 6.0 feet. In figure 10 there is 
no evidence of the beginning of similarity downstream of 
the 0.25-inch rod, and comparison of the profile at 2=10.5 
feet with other distributions, at the same Reynolds number, 
which do exhibit similarity indicated that, even at 2=10.5 
feet, 6.5 feet downstream from the rod, the velocity distri- 
bution is not yet similar. Figures 11 (a) and 11 (b) show 
the distributions for the 0.04-, 0.1 1-, 0.18-, and 0.25-inch 
rods at 2 = 6.0 feet and 2=10.5 feet, respectively. It is seen 
that the deviation among the profiles is less at 2=10.5 feet 
than at 2 = 6.0 feet, illustrating the trend toward similarity 
and the dependence of the similarity position on the rod 
diameter. The trend in figure 10 is in the opposite direction 
of a Reynolds number effect, but the Reynolds number 
change involved is small enough to be negligible compared 
with the differences indicated. It is expected that immedi- 
ately downstream from the rod there would be an influence 
on dP/dx and dP/dy, but there is no evidence that it per- 
sisted to any great length downstream. If the similarity 
position for the 0.04-inch rod is taken at 2=5.5 feet, that is, 
1.5 feet downstream from the rod, the ratio of similarity 
distance to rod diameter is 450. Assuming a linear effect 
of the rod diameter on the similarity distance would put 
the similarity position for the 0.1 1-, 0.1 8-, and 0.25-inch 
rods at 2 = 8 . 2 , 10.9, and 13.5 feet, respectively. This is 
applied only for the condition of the rod at 4 feet and Ui= 108 
feet per second. The wind speed and the position of the 
rod in the layer, that is, the size of the rod compared with 
the bound ary -layer thickness, may be other variables which 
have an effect. The conclusion from this is that too much 
length is required for distortions to disappear for rods to be 
practical thickening devices. 

The velocity distributions at a free-stream speed £A=43 
feet per second for a 24-mesh screen making an angle of 
0.5° with the plate are shown in figure 12. Although no 
summary plot is given, it may be stated that the profiles 
from 2 = 3.5 feet to 2=10.5 feet exhibited similarity. The 
boundary-layer thickness at 2=10.5 feet is approximately 
3 inches as compared with 1.4 inches for the condition of 
free transition. The effect of wider angle or greater solidity 
was to move the similarity position farther downstream and 
to give a thicker layer. For example, the profiles produced 
by a 24-mesh screen, at an angle of 2.5° with the plate, did 
not show similarity at 2=5.0 feet. At 2=9.5 feet the flow 
was similar, and the layer was 5.2 inches thick. Since no 
measurements were made between 2=5.0 feet and 2=9.5 
feet, the position at which similarity first occurred is unknown. 
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Figure 6— Nondimensional mean-velocity distributions for free transition as shown in 
figure 2. U\ = 108 feet per second. 

At an angle of 5° the flow was no longer similar for the length 
of the plate. The use of screens seemed to be an effective 
method of thickening but was not further investigated 
because of the difficulties entailed in mounting a screen, 
keeping it under tension, and maintaining its shape for long 
periods of time. 

Figures 13, 14, and 15 show the velocity distributions 
downstream from the distributed roughness at free-stream 
speeds of 35, 55, and 108 feet per second. Since the general 
characteristics for the three speeds are similar, the nondi- 
mensional profiles are given in figure 16 only for fTi = 55 feet 
per second. It is seen that a distance of about 3.0 feet was 
required for similarity and that at z=10.5 feet the layer was 
approximately 3 inches thick, which for U\= 55 feet per 
second is the thickness of a layer having its beginning about 
14 feet upstream. In other words, the virtual origin of the 
layer was displaced 3.5 feet upstream of the leading edge. 
There was some gain in effective distance at the higher 
speed. At C/i = 108 feet per second, the virtual origin was 
about 5.5 feet upstream of the leading edge. 

Comparison of the velocity distributions in the various 
artificially thickened layers is made with that of free transi- 
tion, in figure 17, at approximately the same Reynolds 
number Re. The agreement is excellent. To test the 
universal character of the velocity distributions, it is of 
interest to compare them with those of other investigations. 
In figure 18, the distributions in the layer thickened by sand 
roughness are compared with the measurements of Hama 
(reference 11) and Schultz-Grunow (reference 12). The 
two curves show the effect of Reynolds number, and the 
data are sufficiently consistent to conclude that the velocity 
distributions obtained for the artificially thickened layer are 


characteristic of a turbulent boundary layer with zero 

pressure gradient and a smooth surface. 


Much of the theoretical work that has been done on the 
turbulent boundary layer is based on the application of the 
logarithmic laws of velocity distribution derived for flow in 
pipes and channels (reference 13). These involve the adop- 
tion of the mixing-length concept, a form for the shearing- 
stress distribution, and the omission of the influence of 
viscosity on the turbulence. 

The wall-proximity law and the velocity-defect law for 
the velocity distribution are compared with experimental 
results in figures 19 and 20. The value of U T in the present 
measurements was determined from the measured velocity 
profile by means of the Von Kdrm&n momentum equation 
which for zero pressure gradient gives 

dd T 0 

dx 2gi 

The mean-velocity distributions for the artificially thickened 
turbulent boundary layer with zero pressure gradient, to- 
gether with those of Schultz-Grunow (reference 12) and 
Freeman (reference 14), are plotted in the form of U/U T 
against yU T lv in figure 19. In contrast with the measure- 
ments for pipes and channels, it is seen that there is a sys- 
tematic increase in the deviation from the logarithmic law as 
the Reynolds number decreases. In the region where the 
law is valid figure 19 gives values of 4.2 for A and 5.75 for 
B. The following table gives the comparison with the 
measurements in smooth pipes and channels for wall 
proximity. 


Investigator 

Type 

A 

B 

Nikuradse (reference 15) - 

Pipe 

5.0 

5.0 

5.0 

7.2 

4.2 

5. 75 
5. 75 
5. 75 
5. 75 
5. 75 

Nlkuradse (reference 16) - 

Qhannel 

Donch (reference 17) 

Channel 

Laufer (reference 5) 

Channel 

Klebanoff and Diehl 

Boundary-layer 



The commonly used values for the experiments of Nikur- 
adse in smooth pipes for points near the wall is 5.5 for B and 
5.8 for A, but Millikan (reference 18) has shown that values 
of 5.75 and 5.0 are in closer agreement with the measure- 
ments. The broken line in figure 19 shows that the power 
law 



is a good approximation for the higher Reynolds number 
R s = 152,000 throughout most of the layer, but it should be 
remarked that at the lower Reynolds numbers the agreement 
was not so good. The distributions at lower Reynolds 
numbers, on a log-log basis, had a sinuous appearance which 
seemed to decrease with increasing Reynolds numbers. 


COMPARISON WITH LOGARITHMIC LAWS 
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FIGURE 8.— Mean-velocity distributions for O.ZWnch cylindrical rod in contact with surface 4 


feet from leading edge. 


l\ = 108 feet per second. 
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Figure 9. — Nondime.nsional mean-velocity distributions for 0.04-inch cylindrical rod in con- 
tact with surface 4 feet from leading edge. I7i = 108 feet per second. 
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0 

Figure 10.— Nondimensional mean-velocity distributions for 0.25-inch cylindrical rod in 
contact with surface 4 feet from leading edge. Ui = 108 feet per second. 
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(a) x=6.0 feet. 

(b) £=10.5 feet. 

Figure 11. —Nondimensional mean-velocity distributions for various rod diameters. L 7 i = 108 

feet per second. 
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Figure 12-Mean-velocity distributions for 1-foot span of 24-mesh screen making angle of 0.5° with plate. U, - 43 feet per second. 
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X 

(ft) 

o 

(in.) 

a 

(in.) 


208 

0.258 

0.145 

7,990 

2.50 

.240 

.157 

8,384 

3.00 

241 

163 

8,720 

400 

.255 

.172 

9,390 

5.00 

.266 

.188 

10,590 

5.50 

.270 

.192 

10,300 

6.50 

294 

.208 

1 1,840 

7.50 

308 

.222 

12,190 

8.50 

.327 

238 

13,060 

950 

.353 

.257 

14,020 

10.50 

.368 

272 

14,850 


1.2 


1.6 

y, in. 


2.0 


2.4 


28 


3.; 


Figure 15.— Mean- velocity distribution? for No. 16 sandpaper covering first 2 feet of surface. Ui = 10S feet per second. 
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Figure 17 —Comparison of nondimensional mean-velocity distributions having approxi- 
mately same values of Ft. 



Figure 16.— Nondimensional mean-velocity distribution for No. 16 sandpaper covering first 
2 feet of surface. Ui = 55 feet per second. 


Figure 18.— Comparison of nondimensional mean- velocity distributions obtained 1 various 

investigators. 
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V 


Figure 19.— Comparison of experimental results with H-power law and logarithmic law. 



Figure 20— Channel and boundary-layer profiles compared on basis of velocity-defect law. 


In figure 20, both boundary-layer and channel results are 
plotted in the form suggested by the velocity-defect law. It 
is seen that results for both the boundary layer and the 
channel fit the straight lines up to y/8= 0.15. The value of 
B seems to be universal for boundary layers and channels 
within this region. While results for pipes are omitted, it is 
remarked that the same holds true for pipe flow. It will be 
noted that the observed results are independent of Reynolds 
number when plotted in this way, indicating separate uni- 
versal relationships of the form 


cient distance downstream, the mean-velocity distributions 
assume a universal character. However, the lack of knowl- 
edge concerning the relationship between the turbulence and 
the mean flow prohibits the use of mean-velocity distribu- 
tions as direct evidence that the turbulence is fully developed. 
Since it is impractical to measure all the characteristics of the 
turbulence and in view of the interest in measurements of the 
spectrum in shear flow, the one-dimensional energy spectrum 
was chosen as the basis of comparison for the turbulence in 
the artificially thickened layers with that resulting from free 
transition. 

The mean-square value of the ^-fluctuations was measured 
between frequencies n and n+dn at values of n from 10 
to 10,000 cycles per second. Over most of the frequency 
range the accuracy of the spectrum measurements is within 
±10 percent. The accuracy was somewhat less at the low- 
frequency end below 100 cycles per second because of large- 
amplitude fluctuations and at the high end above 6000 
cycles per second because of noise and possible wire-length 
effect. No wire-length corrections were made, and the 
noise level at 10,000 cycles per second was approximately 
50 percent of the total signal. A direct over-all measure- 
ment of the intensity u' /Ui was made during each spectrum 
determination to compare with that calculated from the 
spectrum and to serve as a check on the accuracy. In 
general, the agreement was within 5 percent. The spectrum 

curves are given in the nondimensional form — against 


nL x 


by 


Here F(n) is the normalized spectrum function defined 


X' 


F(n) dn = 1 


and L x is the integral scale defined by 
L x — ^ R x dx 


( 1 ) 


where R x is the longitudinal space correlation. The value 
of L x was obtained from the intercept of the F(n) curve 
at n=0 on the assumption that a space correlation was 
related to a time correlation by x=Ut and using the Fourier 
transform written in the following form: 


Fin) 


_± r 

“*7jo 


„ 2tmx , 
R x cos — jj— dx 


( 2 ) 


From equations (1) and (2), when 7i=0 


Ui-U 

U T 


-fiylb) 


T UF(n ) H . 0 

A 


(3) 


A point to be made with emphasis is that boundary-layer 
profiles are not like char *xel and pipe flow profiles, and it is 
seen that an artificially thickened boundary layer conforms 
to the pattern established in figures 19 and 20 after the 
similarity stage has been reached. 

SPECTRUM AND INTENSITY OF u-FLUCTUATION 

It has been shown in the previous section that, at a suffi- 


The values of L x calculated by equation (3) for the turbu- 
lent boundary layer resulting from free transition and the 
boundary layers artificially thickened with a 0.04-inch rod 
and by roughness as described earlier are given in figures 
21, 22, 23, and 24. It is noted that, except near the wall, 
L x is constant across the layer and is proportional to the 
boundary-layer thickness, having a value of LJb of approxi- 
mately 0.4. It should be pointed out that the exactness of 
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equation (2) is open to question at high turbulence levels, 
since it requires that the fluctuations at a point be the result 
of a pure translation with velocity U of a fixed turbulence 
pattern. 

For the condition of free transition, figure 21 shows a set 
of curves at x =10.5 feet, with ?7i = 108 feet per second, at 
different distances from the surface. It is seen that there is 
a systematic variation across the boundary layer up to 
^//5 = 0.57 . Since L x was used as the reference length, all 
curves are brought together at zero frequency. The trend 
in going toward the surface is for the higher frequencies 
(smaller eddies) to have proportionately more of the energy. 

It may be that near the surface some of the smaller eddies 
can derive energy from the steep velocity gradient which 
exists and that farther away from the plate where the velocity 
gradient is comparatively small local isotropy begins at a 
lower wave number. It is hoped that future measurements 
of the spectrum of shearing stress will throw some light on 
this point. 

Comparison is made in figure 22, at x= 10.5 feet and 
y/8= 0.03 and y/ 5=0. ‘57, with the spectra obtained in the 
boundary layer resulting from free transition and those 
obtained in the artificially thickened layers using a 0.04- 
inch rod and sand roughness. The free-stream speed was 108 
feet per second for the rod and free transition and 55 feet 
per second for the roughness so that the boundary-layer 
Reynolds number is approximately the same. The spectra 
are similar to one another and the same family of spectra 
exists in the artificially thickened layers as for that resulting 
from free transition. 

Spectrum measurements made at various distances down- 
stream from the 0.04-inch rod, at y/8 values of 0.57 and 0.4 
and C7i = 108 feet per second, are given in figure 23. It will 
be noted that the curves for x=6.0, 7.5, and 10.5 feet are 
similar. Although a single curve is drawn through the 
measured points, there is evidence of a Reynolds number 
effect in the downstream direction on the high-frequency 
end of the spectrum. The energy associated with the higher 
frequencies is so small that it does not distort the normalized 
curve at the low-frequency end. The curve for x=4.17 feet 
at y/8=0A which is only 2 inches from the rod is definitely 
dissimilar. This curve is practically flat up to 1500 cycles 
per second. The same characteristic flatness was also ob- 
served at y/8= 0.03, so no particular significance should be 
attached to the value of y/8=0A, as compared with the rest 
of the curves with y/8= 0.57. Measurement of the spectrum 
for the condition of free transition, at x=6.9 feet, which was 
as close as it was possible to approach the transition because 
of the back-and-forth movement of transition, showed no 
evidence of the characteristic flatness observed close to the 
rod. In fact, it gave evidence that similarity of the spectra 
began at the downstream extreme of the transition region. 
In the case of free transition, the spectrum was also measured 
in the laminar layer at x=4.5 feet. A very sharp peak was 
found at 280 cycles per second, the predominant frequency 
of the laminar-boundary-layer oscillations, and little energy 
was found at other frequencies. There was no evidence of 


this peak in the turbulent boundary layer at z=6.9 feet. 
It was thought at first that the flatness of the spectrum ob- 
served near the rod was associated with the distorted mean- 
velocity profile and that here was evidence of a connection 
between turbulence distortion and mean-velocity distortion. 
However, figure 24, which pertains to thickening by rough- 
ness, shows no effect in the spectrum at the 3-foot position 
where the mean-velocity profile was definitely distorted. It 
would seem that the flatness of the spectrum measured near 
the rod is an* effect peculiar to the rod and may be associated 
with eddy frequencies from the rod superimposed on a more 
or less normal spectrum. 

It may be expected that, in the turbulent boundary layer, 
at a distance far enough away from the wall an extensive 
region of local isotropy exists. It is of interest to see how 
well the spectra in the artificially thickened turbulent 
boundary layer where the Reynolds number of turbulence 
is high agree with existing theories predicting the shape of 
the spectrum in isotropic turbulence. 

Before making the comparison it may be well to discuss 
briefly some important features of the theory. A funda- 
mental concept in the theories introduced by Kolmogoroff 
(references 1 and 2) and Heisenberg (reference 19) is that 
energy passes into the spectrum by way of the larger eddies 
and is then passed on to smaller and smaller eddies. Ap- 
preciable dissipation by viscosity starts only when the eddies 
get small enough. Beyond a certain wave number, smaller 
eddies get their energy only from larger eddies and are sup- 
plied sufficient energy to be kept in statistical equilibrium 
even in cases where the larger eddies are decaying. When 
the Reynolds number of turbulence u'L x lv is high, there will 
exist in the lower-wave-number part of the equilibrium 
range a subrange in which inertial forces predominate and 
the transfer of energy is equal to the total dissipation e and 
at the high-wave-number end another subrange where viscous 
effects predominate. Heisenberg assumed that the transfer 
of energy at wave number k was caused bv ? a turbulence 
friction produced by eddies with wave numbers greater 
than k and that the total dissipation involved a turbulence 
viscosity plus ordinary viscosity in an expression of the 
following form: 

t= [" +x lVw dk ^l^ k y E s k, ) dk ' ( 4 ) 

Here the second term within the brackets represents the 
turbulence viscosity, and x is an absolute constant. The 
solution of equation (4) gives the spectrum in the equilibrium 
range the form 

*<K0' V ‘1 1+ (£)T'’ (5) 

where k s is a wave number in the intermediate range where 
inertial forces and viscosity both play a role and is given by 
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Figure 22.— Spectrum of u-fluctuations in turbulent boundary layer at .r= 10.5 feet. U\ = 108 feet per second for free transition and 0.04-inch rod; £71 = 55 feet per second for roughness. 
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Figure 24.— Test of -5/3 law and -7 law for spectrum 


of u-fluctuations in turbulent boundary layer thickened by roughness for first 2 feet of surface. 


(7, =55 feet per second; j//5=0.57. 
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The equilibrium spectrum given by equation (5) deals 
with a three-dimensional energy spectrum. In order to 
compare with experiment it is necessary to transform this 
equation to the one-dimensional form using the following 
relationship also derived by Heisenberg: 




The one-dimensional form of equation (5) is found to be 


E x (k x )-- 


110 


Aik x 


where 


and 


■N&r 


( 6 ) 


(7) 


Ai 


/8A 2/3 

\9x/ 


Ic sx =0.645k s 

In the inertial range, that is, equation (7) becomes 


^ ( ** )= TTo ^*~ 5/3 (8) 

and in the viscous range, that is, k x ^k s , equation (7) becomes 


where 


E x {k x ) 


1 

126 


A 2 k x ~ 7 



(9) 


In figure 24, equations (8) and (9) are compared with the 
spectra measured at different distances downstream from 
the roughness, at y/8= 0.57 and £/i = 55 feet per second. At 
the high-turbulence Reynolds number of the boundary layer, 
the existence of a —5/3 region is quite evident. The agree- 
ment with the n~ 7 law should be taken with some reservation 
because of the lesser accuracy of the measurements at the 
high-frequency end of the spectrum. It should be men- 
tioned in passing that the slope of the spectrum curves at 
y/8= 0.03, such as shown in figure 22, is approximately 1, 
Where the —5/3 region exists at y/8= 0.57. 

It is seen from figure 25 that in the intermediate range of 
frequencies equation (7) does not agree so well as the inter- 
polation formula 


£.(fc)- I f 5 AI:,-'-[l+(^)*' , ]- ! (10) 

suggested by Heisenberg. For convenience, equations (7) 
and (10) are represented in figure 25 as 


and 


where 


and 


^>--©‘T+(£)T 


F(z) 


UF(n) 

L x 


(7a) 

(10a) 


and A 0 and z a are experimentally determined. Formula 
(10a), as shown in figure 26, is also in good agreement with 
some unpublished measurements made in isotropic turbulence 
50 inches downstream from a 1-inch grid in the NBS 4^-foot 
wind tunnel. The fact that the analytical relation (7) does 
not agree with the measured spectra so well as the empirical 
expression (10) would indicate that equation (4) may not yet 
be quite exact. Expressions (7a) and (10a) are identical in 
the extreme regions of the equilibrium range, and it is of 
interest to calculate the value of x from the constants A 0 and 
2 *. Since A 0 and 2 , are determined from a nondimensional 
normalized representation they are connected to x by the 
following relations. 

The term E{k) of the three-dimensional spectrum is defined 
such that 

rE(k)dk J ?± ?+'* 

Jo 2 

From relation (6), it can be shown that 

j o °°E x (k x ) dk x =^CE (it) dk 
so that, for isotropic turbulence where u 2 = v 2 =w 2 , 


and 


since 


j; 


E x (k x )dk x ,=4 

0 4 


F(z)dz = J“E x (k x )dk x 


2i m 27 vz 




In the inertial range, 


and 


F(z)=A 0 z s 5 '*z ~ 5/3 

1 /27r\ 5/3 

E x (k x ) = A 0 z, m L x u 2 \j^) k x 


- 5/3 




5/3 


(ii) 


It will be noted that e and u 2 must be determined in 
addition to the quantities already mentioned before x can 
be calculated by equation (11). For isotropic turbulence, 
the expression for the dissipation (reference 20) is 

e=15 "^ ( 12 ) 


where X is the microscale obtainable from a spectrum 
measurement by 


nL : 

Z —U 


5?-£f *'*<">*• 


(13) 


ARTIFICIALLY THICKENED FULLY DEVELOPED TURBULENT BOUNDARY LAYERS WITH ZERO PRESSURE GRADIENT 23 



24 


REPORT 1 1 10 — 1 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



FI 0VKK 26. — Spectrum »r isotropic turbuience compared with Ration ConnuJ* 0 “ " ^ '' inCh ^ ^ ““ ^ 

u'/U= 0.013; £7=50 feet per second; L, =0.047 foot,_and > -o.oio loot 


= 0 . 0162 ; 


z.= 4.00; 
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Partial justification for the use of such relationships where 
the turbulence does not have general isotropy is afforded by 
figure 27, where it is seen that about 90 percent of the integral 
in equation (13) is contributed by frequencies above 500 
cycles per second. This is in the region where the —5/3- and 
— 7- power laws hold and is th eref ore the region of local isotropy. 
It is concluded therefore as a reasonable approximation that 
X determined by equation (13) is the appropriate and only 
scale entering into the viscous dissipation. However, much 
of the value of u 2 is contributed. below 500 cycles per second, 
and the assumption that u 2 -{-v 2J rw 2 equals 3 u 2 , as made in 
equation (12), is not a sufficiently good approximation. 
Since v 2 and vP were not measured, the results of reference 6 
were used to obtain 


u 2 +v 2 +w 2 =0.62(3u 2 ) 

at y/8= 0.57. The factor 0.62 is applied to_the right side of 
equation (12) and to equation (11) where u 2 appears. 

The value of x calculated from equation (11) for the 
turbulent boundary layer artificially thickened by sand 


roughness, at 2= 10.5 feet and y/8= 0.57, with the values of 
A 0 and z s given in figure 25, is 0.23. The value of x calcu- 
lated for the case of isotropic turbulence shown in figure 26 
is 0.65. By a similar procedure x was calculated for the 
spectrum measured by Laufer in fully developed channel 
flow (reference 5), and a value of 0.52 was obtained. At 
present there are insufficient experimental data to determine 
the reason for the variations obtained. It is of interest 

u'L 

to mention that the turbulence Reynolds numbers were 


1790, 550, and 190 for the boundary layer, channel, and 
isotropic turbulence, respectively, and that corresponding 
increasing values of x were obtained. 

The insensitivity of the spectrum to distortions in the 
mean- velocity distributions as illustrated by figure 24 
made it desirable to check further on the fully developed 
nature of the turbulence in the artificially thickened layer 
by measuring the distribution of u'/Ui. Intensity measure- 
ments were made at various distances downstream from the 
roughness at C7i = 108 feet per second, and these are given 
in figure 28. At z=3.0 feet, 1 foot downstream from the 


rtF{n) 
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edge of the roughness, the intensity is flat up to y 1 8=0 A 
and is apparently the result of the diffusion of the high 
turbulence generated by the roughness. From x=5.0 to 
2=10.5 feet, the distributions of u' exhibit similarity, which 
is the same region for which the mean-velocity distributions 
were similar. The change in turbulence level accompanying 
a change in speed from 108 to 55 feet per second is also shown. 
It should be noted that there is no change in level accom- 
panying the variations in boundary-layer Reynolds number 
due to changes in boundary-layer thickness. 

An interesting feature of the turbulence in the boundary 
layer was its intermittent nature in the outer portions of the 
layer. This was also observed by Townsend for a wake 
(reference 21) and by Corrsin for a turbulent jet (reference 
22). It seems to be characteristic of shear flow with a free 
boundary. Figure 29 shows two oscillograms taken at 
z=10.5 feet, at ?//5=0.03 and 2//^= 0.75. At y/8= 0.03 there 



Figure 28— Intensity of u-component of turbulence for boundary layer thickened by No. 16 
sandpaper over first 2 feet of surface. 


is no trace of intermittency, while at y/8= 0.75 it is quite 
marked. Observation of the oscilloscope trace indicated 
there is a gradual increase in the intermittency as the distance 
from the surface increased. It is planned in the future to 
make some quantitative measurements of the intermittency 
and its effect on the measured turbulence intensity 

CONCLUSIONS 

From tests made to determine the feasibility of artificially 
thickening a turbulent boundary layer on a flat plate, the 
following conclusions can be drawn: 

1. It has been established that by artificial methods it is 
possible to thicken substantially a fully developed turbulent 
boundary layer on a flat plate. The layer so formed be- 
comes characteristic of that developed from the beginning 
on a smooth flat plate after passing over the required length 
of smooth surface. For a given boundary -layer thickness, 
artificial thickening makes possible a saving in length pro- 
vided that there is enough length for the layer to become 
free from distortions. 

2. Almost any kind of obstruction placed on the surface 
would produce a thickened turbulent flow, but only from a 
cerfain class did downstream effects disappear in a reason- 
ably short distance. This class appeared to be character- 
ized by the ability to absorb energy from the mean flow 
without introducing large-scale disturbance. Two examples 
of this class are as follows: 

(1) A fine-mesh screen placed so that the plate and screen 

form a narrow wedge opening into the wind 

(2) Sand roughness cemented to the surface 

3. The mean- velocity distribution appears to afford a 
simple means of identifying a fully developed turbulent 
boundary layer. 

4. The logarithmic laws of velocity distribution are not in 
good agreement with the measured distributions in the 
turbulent boundary layer with zero pressure gradient and 
are applicable only in a small region near the wall. The 
velocity profiles in the boundary layer are independent of 
Reynolds number on a velocity-defect basis, indicating a 


A/#\J 

(O) 



(a) y/S= 0.03. 

(b) y IS =0.75. 

Figure 20. Oscillograms of w-fluctuations in turbulent boundary layer. t= 10.5 feet;Ui=55 feet per second; 5=3.1 inches. Time interval between dots, Ho second. Boundary layer thick- 
ened by No. 16 sandpaper over first 2 feet of surface. 
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universal relationship of the form 


where U is the mean velocity in the boundary layer, U\ is 
the mean velocity in the free stream just outside the bound- 
ary layer, U T is the square root of the ratio of the shearing 
stress at the surface to the density, y is the distance normal 
to the surface measured from, the surface of a flat plate, and 
5 is the boundary-layer thickness. 

5. The longitudinal scale of turbulence L x in a fully de- 
veloped layer is proportional to the boundary-layer thickness 
5, with a value of L x /b of approximately 0.4, and except near 
the wall is constant across the layer. 

6. There is a systematic variation in the one-dimensional 
energy spectrum across the boundary layer up to y/b= 0.6, 
and in the outer region of the boundary layer, that is, 
i//5> 0.6, the existing theories for the spectrum in isotropic 
turbulence can be applied. 

7. There is an intermittent character to the flow in the 
turbulent boundary layer similar to that found for wakes 
and jets. 

National Bureau of Standards, 

Washington, D. C., September 15, 1950 . 
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